Um estudo da atenção primária à saúde mediante o indicador: Internações por condições sensíveis à atenção ambulatorial (CSAA)

A Atenção Primária à Saúde - APS é reconhecida como o nível fundamental e porta de entrada do sistema de atenção à saúde, sendo o lugar adequado onde pode ser atendida e resolvida a maior parte dos problemas de saúde. É considerada pela OMS como a principal proposta de modelo assistencial. Essa importância da APS leva a necessidade de pesquisas avaliativas dos seus resultados para adequação e melhoria de políticas e planos de ação delineados em relação à mesma. Pesquisas internacionais e nacionais são realizadas, nas quais indicadores relativos às atividades hospitalares estão sendo empregados com o objetivo de medir resultados como efetividade e acesso da APS. Um desses indicadores, desenvolvido por John Billings da Universidade de Nova York, na década de 90, consiste nas condições pelas quais as internações hospitalares por Condições Sensíveis à Atenção Ambulatorial (CSAA) deveriam ser evitadas caso os serviços da APS fossem efetivos e acessíveis. Utilizando-se o SIH-AIH/2008 e a lista brasileira de Internações por Condições Sensíveis a Atenção Primária, publicada em 2008, a proposta do presente trabalho é a de estudar os cuidados primários à saúde baseando-se nas ICSAA, na área urbana da cidade de Juiz de Fora-MG. Buscou-se responder sobre os efeitos que ocorrem nessas internações a partir das características individuais dos pacientes, das características das Unidades Básicas de Saúde-UBS (infraestrutura, produção e modelos assistenciais) e das condições sócio-econômicas/ambientais das áreas cobertas por UAPS e descobertas (sem UAPS), com a utilização de modelos multiníveis logísticos com intercepto aleatório. Buscou-se conhecer, também, a distribuição espacial das taxas padronizadas por idade das ICSAA nessas áreas e suas associações com as variáveis contextuais, utilizando-se ferramentas da análise espacial. Os resultados do presente trabalho mostraram que a porcentagem de internações por CSAA, foi de 4,1%. Os modelos assistenciais ESF e o Modelo Tradicional, base da organização da atenção primária no Brasil, não apresentaram no município, impacto significativo nas ICSAA, somente na forma de áreas descobertas tendo como referência as áreas cobertas. Também não foram significativas as variáveis de infraestrutura e produção das UAPS. Os efeitos individuais (idade e sexo) nas ICSAA foram significativos, apresentando probabilidades de significância menores que 1%, o mesmo acontecendo com o Índice de Desenvolvimento Social-IDS, que contempla as condições sociais, econômicas e ambientais das áreas analisadas. A distribuição espacial das taxas padronizadas por idade apresentou padrão aleatório e os testes dos Multiplicadores de Lagrange não foram significativos indicando o modelo de regressão clássico (MQO) como adequado para explicar as taxas em função das variáveis contextuais. Para a análise conjunta das áreas cobertas e descobertas foram fatores de risco: a variável econômica (% dos domicílios com renda até 2 SM), áreas descobertas tendo como referência as áreas cobertas e a região nordeste do município. Para as áreas cobertas as variáveis de produção das UAPS, econômica e a região nordeste apresentaram como fator de risco para as taxas de internação por CSAA.

Tópicos em programação dinâmica

Muitos problemas de Dinâmica em Economia se encaixam dentro de uma estrutura de modelos de decisão seqüencial, sendo resolvidos recursivamente. Programação Dinâmica uma técnica de otimização condicionada que se encarrega de solucionar problemas desse tipo. Esse trabalho tem como objetivo apresentar uma resenha dos principais resultados teóricos em Programação Dinâmica. Os métodos da Programação Dinâmica são válidos tanto para problemas determinísticos como para os que incorporam variável incerteza. esperada objetividade de uma dissertação de Mestrado, no entanto, nos impediu de extender análise, deixando assim de considerar explicitamente neste trabalho modelos estocásticos, que teria enriquecido bastante parte destinada aplicações Teor ia Econômica. No capítulo desenvolvemos instrumental matemático, introduzindo uma série de conceitos resultados sobre os quais se constrói análise nos capítulos subsequentes. Ilustramos tais conceitos com exemplos que seguem um certo encadeamento. Nas seções 1.1 1.2 apresentamos as idéias propriedades de espaços métricos espaços vetoriais. Na seção 1.3, prosseguimos com tópicos em análise funcional, introduzindo noção de norma de um vetor de espaços de Banach. seção 1.4 entra com idéia de contração, Teor ema do Ponto Fixo de Banach e o teor ema de Blackwell. O Teorema de Hahn-Banach, tanto na sua forma de extensão quanto na sua forma geométrica, preocupação na seção 1.5. Em particular, forma geométrica desse teorema seus corolários são importantes para análise conduzida no terceiro capítulo. Por fim, na seção 6, apresentamos Teorema do Máximo. Ao final deste capítulo, como também dos demais, procuramos sempre citar as fontes consultadas bem como extensões ou tratamentos alternativos ao contido no texto. No capítulo II apresentamos os resultados métodos da Programação Dinâmica em si seção 2.1 cuida da base da teoria, com Princípio da Otimal idade de Eellman e a derivação de um algoritmo de Programação Dinâmica. Na seção 2.2 mostramos que esse algoritmo converge para função valor ótima de um problema de horizonte infinito, sendo que esta última satisfaz chamada Equação de Bellman. seção seguinte se preocupa em fornecer caracterizaçBes para função valor mencionada acima, mostrando-se propriedades acerca de sua monotonicidade concavidade. seção 2.4 trata da questão da diferenciabi idade da função valor, que permite se obter alguns resultados de estática Cou dinâmica} comparativa partir da Equação de Bellman. Finalmente, na seção 2.5 apresentamos uma primeira aplicação Teoria Econômica, através de um modelo de crescimento econômico ótimo. No capítulo III introduzimos uma outra técnica de otimização Programação Convexa- mostramos dificuldade em se tentar estabelecer alguma relação de dominância entre Programação Dinâmica Programação Convexa. Na seção 3.2 "apresentamos os Teoremas de Separação, dos quais nos utilizamos na seção seguinte para demonstrar existência de Multiplicadores de Lagrange no problema geral da Programação Convexa. No final desta seção dizemos porque não podemos inferir que em espaços de dimensão infinita Programação Convexa não pode ser aplicada, ao contrário da Programação Dinâmica, que evidenciaria uma dominancia dessa última técnica nesses espaços. Finalmente, capítulo IV destinado uma aplicação imediata das técnicas desenvolvidas principalmente no segundo capítulo. Com auxílio dessas técnicas resolve-se um problema de maximização intertemporal, faz-se uma comparação dos resultados obtidos através de uma solução cooperativa de uma solução não-cooperativa.

Um método primal-dual aplicado na resolução do problema de fluxo de potência ótimo

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Análise não suave e aplicações em otimização

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Convexidade generalizada em problemas de controle ótimo com tempo livre

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Economic optimization method to design telescope irrigation of multiples outlets

In this study is presented an economic optimization method to design telescope irrigation laterals (multidiameter) with regular spaced outlets. The proposed analytical hydraulic solution was validated by means of a pipeline composed of three different diameters. The minimum acquisition cost of the telescope pipeline was determined by an ideal arrangement of lengths and respective diameters for each one of the three segments. The mathematical optimization method based on the Lagrange multipliers provides a strategy for finding the maximum or minimum of a function subject to certain constraints. In this case, the objective function describes the acquisition cost of pipes, and the constraints are determined from hydraulic parameters as length of irrigation laterals and total head loss permitted. The developed analytical solution provides the ideal combination of each pipe segment length and respective diameter, resulting in a decreased of the acquisition cost.

Comparative Analysis of Tolerances Allocation in Mechanical Assemblies based on Cost-tolerance Curves

El presente trabajo propone un método para la determinación de los valores de las tolerancias individuales de las piezas que forman un conjunto ensamblado a partir de valores de tolerancias especificados en el conjunto final, optimizando el coste total de fabricación de las piezas individuales a partir de funciones de coste-tolerancia basadas en el proceso de fabricación de cada una de ellas. Para ello se parte de los principales trabajos desarrollados en la línea de asignación de tolerancias y se realiza la propuesta del modelo de trabajo, basado en la optimización de costes a partir de la aplicación del método de los multiplicadores de Lagrange a diversas curvas de coste-tolerancia

Tratamiento consistente de restricciones e impacto en dinámica de sistemas multicuerpo

Esta tesis está enmarcada en el estudio de diferentes procedimientos numéricos para resolver la dinámica de un sistema multicuerpo sometido a restricciones e impacto, que puede estar compuesto por sólidos rígidos y deformables conectados entre sí por diversos tipos de uniones. Dentro de los métodos numéricos analizados se presta un especial interés a los métodos consistentes, los cuales tienen por objetivo que la energía calculada en cada paso de tiempo, para un sistema mecánico, tenga una evolución coherente con el comportamiento teórico de la energía. En otras palabras, un método consistente mantiene constante la energía total en un problema conservativo, y en presencia de fuerzas disipativas proporciona un decremento positivo de la energía total. En esta línea se desarrolla un algoritmo numérico consistente con la energía total para resolver las ecuaciones de la dinámica de un sistema multicuerpo. Como parte de este algoritmo se formulan energéticamente consistentes las restricciones y el contacto empleando multiplicadores de Lagrange, penalización y Lagrange aumentado. Se propone también un método para el contacto con sólidos rígidos representados mediante superficies implícitas, basado en una restricción regularizada que se adaptada adecuadamente para el cumplimiento exacto de la restricción de contacto y para ser consistente con la conservación de la energía total. En este contexto se estudian dos enfoques: uno para el contacto elástico puro (sin deformación) formulado con penalización y Lagrange aumentado; y otro basado en un modelo constitutivo para el contacto con penetración. En el segundo enfoque se usa un potencial de penalización que, en ausencia de componentes disipativas, restaura la energía almacenada en el contacto y disipa energía de forma consistente con el modelo continuo cuando las componentes de amortiguamiento y fricción son consideradas. This thesis focuses on the study of several numerical procedures used to solve the dynamics of a multibody system subjected to constraints and impact. The system may be composed by rigid and deformable bodies connected by different types of joints. Within this framework, special attention is paid to consistent methods, which preserve the theoretical behavior of the energy at each time step. In other words, a consistent method keeps the total energy constant in a conservative problem, and provides a positive decrease in the total energy when dissipative forces are present. A numerical algorithm has been developed for solving the dynamical equations of multibody systems, which is energetically consistent. Energetic consistency in contacts and constraints is formulated using Lagrange multipliers, penalty and augmented Lagrange methods. A contact methodology is proposed for rigid bodies with a boundary represented by implicit surfaces. The method is based on a suitable regularized constraint formulation, adapted both to fulfill exactly the contact constraint, and to be consistent with the conservation of the total energy. In this context two different approaches are studied: the first applied to pure elastic contact (without deformation), formulated with penalty and augmented Lagrange; and a second one based on a constitutive model for contact with penetration. In this second approach, a penalty potential is used in the constitutive model, that restores the energy stored in the contact when no dissipative effects are present. On the other hand, the energy is dissipated consistently with the continuous model when friction and damping are considered.

Efectos dinámicos debidos al tráfico de ferrocarril sobre la infraestructura de vía y las estructuras

La vía tradicional sobre balasto sigue siendo una selección para las líneas de alta velocidad a pesar de los problemas técnicos y la prestación del funcionamiento. El problema de la vía sobre balasto es el proceso continuo del deterioro de éste debido a las cargas asociadas al tráfico ferroviario. En consecuencia es imprescindible un mantenimiento continuado para mantener un alineamiento adecuado de la vía. Por eso se surge la necesidad de comprender mejor el mecanismo involucrado en el deterioro de la vía y los factores claves que rigen su progresión a lo largo de ciclos de carga con el fin de reducir los costos del mantenimiento de la vía y mejorar el diseño de las nuevas vías. La presente tesis intenta por un lado desarrollar los modelos más adecuados y eficientes del vehículo y de la vía para los cálculos de los efectos dinámicos debido al tráfico de ferrocarril sobre la infraestructura de la vía sobre balasto, y por otro evaluar estos efectos dinámicos sobre el deterioro de la vía sobre balasto a largo plazo, empleando un adecuado modelo de predicción del deterioro de la misma. Se incluye en el trabajo una recopilación del estado del arte en lo referente a la dinámica de la vía, a la modelización del vehículo, de la vía y de la interacción entre ambos. También se hace un repaso al deterioro de la vía y los factores que influyen en su proceso. Para la primera línea de investigación de esta tesis, se han desarrollado los diferentes modelos del vehículo y de la vía y la modelización de la interacción entre ambos para los cálculos dinámicos en dos y tres dimensiones. En la interacción vehículo-vía, se ha empleado la formulación de contacto nodo-superficie para establecer la identificación de las superficies en contacto y el método de los multiplicadores de Lagrange para imponer las restricciones de contacto. El modelo de interacción se ha contrastado con los casos reportados en la literatura. Teniendo en cuenta el contacto no lineal entre rueda-carril y los perfiles de irregularidades distribuidas de la vía, se han evaluado y comparado los efectos dinámicos sobre el sistema vehículo-vía en la interacción de ambos, para distintas velocidades de circulación del vehículo, en los aspectos como la vibración del vehículo, fuerza de contacto, fuerza transmitida en los railpads, la vibración del carril. También se hace un estudio de la influencia de las propiedades de los componentes de la vía en la respuesta dinámica del sistema vehículo-vía. Se ha desarrollado el modelo del asiento de la vía que consiste en la implementación del modelo de acumulación de Bochum y del modelo de hipoplasticidad en la subrutina del usuario \UMAT" del programa ABAQUS. La implementación numérica ha sido comprobado al comparar los resultados de las simulaciones numéricas con los reportados en la literatura. Se ha evaluado la calidad geométrica de la vía sobre balasto de los tramos de estudio con datos reales de la auscultación proporcionados por ADIF (2012). Se ha propuesto una metodología de simulación, empleando el modelo de asiento, para reproducir el deterioro de la geometría de la vía. Se usan los perfiles de la nivelación longitudinal de la auscultación como perfiles de irregularidades iniciales de la vía en las simulaciones numéricas. También se evalúa la influencia de la velocidad de circulación sobre el deterioro de la vía. The traditional ballast track structures are still being used in high speed railways lines with success, however technical problems or performance features have led to ballast track solution in some cases. The considerable maintenance work is needed for ballasted tracks due to the track deterioration. Therefore it is very important to understand the mechanism of track deterioration and to predict the track settlement or track irregularity growth rate in order to reduce track maintenance costs and enable new track structures to be designed. This thesis attempts to develop the most adequate and efficient models for calculation of dynamic track load effects on railways track infrastructure, and to evaluate these dynamic effects on the track settlement, using a track settlement prediction model, which consists of the vehicle/track dynamic model previously selected and a track settlement law. A revision of the state of the knowledge regarding the track dynamics, the modelling of the vehicle, the track and the interaction between them is included. An overview related to the track deterioration and the factors influencing the track settlement is also done. For the first research of this thesis, the different models of vehicle, track and the modelling of the interaction between both have been developed. In the vehicle-track interaction, the node-surface contact formulation to establish the identification of the surfaces in contact and the Lagrange multipliers method to enforce contact constraint are used. The interaction model has been verified by contrast with some benchmarks reported in the literature. Considering the nonlinear contact between wheel-rail and the track irregularities, the dynamic effects on the vehicle-track system have been evaluated and compared, for different speeds of the vehicle, in aspects as vehicle vibration, contact force, force transmitted in railpads, rail vibration. A study of the influence of the properties of the track components on the the dynamic response of the vehicle-track system has been done. The track settlement model is developed that consist of the Bochum accumulation model and the hipoplasticity model in the user subroutine \UMAT" of the program ABAQUS. The numerical implementation has been verified by comparing the numerical results with those reported in the literature. The geometric quality of the ballast track has been evaluated with real data of auscultation provided by ADIF (2012). The simulation methodology has been proposed, using the settlement model for the ballast material, to reproduce the deterioration of the track geometry. The profiles of the longitudinal level of the auscultation is used as initial profiles of the track irregularities in the numerical simulation. The influence of the running speed on the track deterioration is also investigated.

Enfoque al ahorro de energía en la conducción de ferrocarriles utilizando el método múltiple de restricción de la optimización

El ahorro de la energía consumida no es un tema menor en el sector ferroviario. Las infraestructuras ferroviarias y los planes operativos implicados en distintos sistemas como el metropolitano, trenes locales, larga distancia y alta velocidad, deben ser realizados teniendo en cuenta las múltiples alternativas para alcanzar un sistema de optimización en el ahorro de energía. Para conseguir alcanzar la mejor solución final, las alternativas deben recoger y trabajar con el máximo de variables asociadas al movimiento del tren y a los planes de operación. El método presentado aquí usa una solución semianalítica que conlleva una discretización del problema, para posteriormente aplicar el método de los Multiplicadores de Lagrange para resolver la optimización de las n-tuplas de velocidad. Con esta alternativa propuesta, es posible incluir todos los detalles sobre las operaciones del tren, como restricciones en el horario o tipos de frenado.

Controlador integrado baseado em conceito híbrido e controle geométrico para embarcações posicionadas dinamicamente.

O conceito de controle híbrido é aplicado à operação de alívio entre um FPWSO e um navio aliviador. Ambos os navios mantêm suas posições e aproamentos pelo resultado da ação do seu Sistema de Posicionamento Dinâmico (SPD). O alívio dura cerca de 24 horas para ser concluído. Durante este período, o estado de mar pode se alterar e os calados estão sendo constantemente alterados. Um controlador híbrido é projetado para permitir modificacões dos parâmetros de controle/observação se alguma alteração significante do estado de mar e/ou calado das embarcações ocorrer. O principal objetivo dos controladores é manter o posicionamento relativo entre os navios com o intuito de evitar perigosa proximidade ou excesso de tensão no cabo. Com isto em mente, uma nova estratégia de controle que atue integradamente em ambos os navios é proposta baseda em geometria diferencial. Observadores não lineares baseados em passividade são aplicados para estimar a posição, a velocidade e as forças externas de mares calmos até extremos. O critério para troca do controle/observação é baseado na variação do calado e no estado de mar. O calado é assumido conhecido e o estado de mar é estimado pela frequência de pico do espectro do movimento de primeira ordem dos navios. Um modelo de perturbação é proposto para encontrar o número de controladores do sistema híbrido. A equivalência entre o controle geométrico e o controlador baseado em Multiplicadores de Lagrange é demonstrada. Assumindo algumas hipóteses, a equivalência entre os controladores geométrico e o PD é também apresentada. O desempenho da nova estratégia é avaliada por meio de simulações numéricas e comparada a um controlador PD. Os resultados apresentam muito bom desempenho em função do objetivo proposto. A comparação entre a abordagem geométrica e o controlador PD aponta um desempenho muito parecido entre eles.

Optimización (curso 2011-2012)

Material docente de la asignatura «Simulación y Optimización de procesos químicos». Parte de Optimización OPTIMIZACIÓN TEMA 6. Conceptos Básicos 6.1 Introducción. Desarrollo histórico de la optimización de procesos. 6.2 Funciones y regiones cóncavas y convexas. 6.3 Optimización sin restricciones. 6.4 Optimización con restricciones de igualdad y desigualdad. Condiciones de optimalidad de Karush Khun Tucker 6.5 Interpretación de los Multiplicadores de Lagrange. TEMA 7. Programación lineal 7.1 Introducción. Planteamiento del problema en forma canónica y forma estándar. 7.2 Teoremas de la programación lineal 7.3 Resolución gráfica 7.4 Resolución en forma de tabla. El método simplex. 7.5 Variables artificiales. Método de la Gran M y método de las dos fases. 7.6 Conceptos básicos de dualidad. TEMA 8. Programación no lineal 8.1 Repaso de métodos numéricos de optimización sin restricciones 8.2 Optimización con restricciones. Fundamento de los métodos de programación cuadrática sucesiva y de gradiente reducido. TEMA 9. Introducción a la programación lineal y no lineal con variables discretas. 9.1 Conceptos básicos para la resolución de problemas lineales con variables discretas.(MILP, mixed integer linear programming) 9.2 Introducción a la programación no lineal con variables continuas y discretas (MINLP mixed integer non linear programming) 9.3 Modelado de problemas con variables binarias: 9.3.1 Conceptos básicos de álgebra de Boole 9.3.2 Transformación de expresiones lógicas a expresiones algebraicas 9.3.3 Modelado con variables discretas y continuas. Formulación de envolvente convexa y de la gran M.

Funciones de varias variables y sus derivadas. Con aplicaciones a la administración, a las Ciencias Sociales y a la Economía

Funciones de varias variables y sus derivadas con aplicaciones a la administración, las Ciencias Sociales y la Economía es una lección en la que, además de los conceptos fundamentales fe funciones de varias variables, se construyen gráficas en dos y tres dimensiones, algunas de ellas utilizando software matemático; de igual manera, se muestran aplicaciones de funciones de producción, costo y demanda, curvas de indiferencia e isocuantas, aplicaciones a la administración, las ciencias sociales y la economía de la diferencia total, análisis marginal, costo promedio marginal, tasa marginal de sustitución, artículos sustitutos y complementarios, y finalmente, se le da importancia al planteamiento de problemas de optimización empleando los multiplicadores de Lagrange.

Optimización de una Función de Utilidad Para un Hogar Promedio de El Salvador.

En esta investigación se esboza una breve historia de cómo y porque surge la teoría de la elección racional, además de exponer los principios generales aplicables a las Ciencias que estudian el comportamiento de la interacción humana, es decir, las Ciencias Sociales. Dichos principios se recogen en economía bajo la rama conocida como microeconomía, específicamente en la teoría de la conducta del consumidor. Es por ello que esta teoría es el pilar fundamental de esta investigación. Dicha teoría se retoma para poder aplicarla al hogar salvadoreño promedio visto como la unidad de consumo dentro de la economía, ya que con la teoría de la conducta del consumidor es posible medir su bienestar (conocida en microeconomía como utilidad) a través la variable proxy consumo mediante la estimación de una función de consumo. Para ello se realizó un diagnóstico del consumo de los hogares salvadoreños en el que se determinó la importancia que tienen en la dinámica económica del país, resultando ser el que mayor peso tiene respecto a los diferentes componentes que conforman el Producto Interno Bruto (PIB). Para tal propósito se realizan tres ejercicios: El primero es una comparación gráfica de cada componente con el PIB a nivel de tasas de crecimiento, en el segundo ejercicio se determinó el porcentaje de contribución de cada componente en el crecimiento del PIB, mientras que en el tercero se realiza una regresión simple para conocer el grado de correlación que tiene cada componente con el PIB. Dicho diagnóstico también permitió conocer la estructura de consumo de los hogares, y cómo esta se ve afectada por los niveles de ingreso e inflación, se determina en grado de importancia de mayor a menor la siguiente estructura de cestas de bienes y consumo constituido por: Alimentos, servicios de vivienda, educación, bienes y servicios diversos, y, salud. Se construyeron las gráficas de Engels y Demanda-Precio para cada canasta de bienes y servicios, con el objetivo de conocer el comportamiento que ha tenido el consumidor ante los precios y los ingresos. Una vez conocida la estructura de consumo se procedió a estimar la función de utilidad del hogar salvadoreño, mediante el método de regresión lineal múltiple, la cual toma la forma de una función del tipo Cobb-Douglas. Estimada la función de utilidad el siguiente que se realizó fue optimizar los niveles de consumo, por el método de los multiplicadores de Lagrange, lo que permitió realizar comparaciones entre los niveles observados de consumo y los niveles óptimos que el hogar puede haber realizado de acuerdo a sus niveles de consumo, lo que revela que los hogares se encuentran consumiendo por debajo de sus niveles óptimos. Bajo un escenario en donde los hogares logren optimizar su consumo, el nivel de bienestar más alto se encuentra en el año 2000 y para alcanzar dicho nivel de utilidad en el año 2013 los hogares necesitan tener ingresos de $672 aproximadamente. Lo anterior se realizó bajo un contexto sin impuestos de ningún tipo, sin embargo, al considerar el impuesto sobre la renta, el óptimo del consumidor sufre un impacto negativo entre el 3% y 5% En base a lo expuesto se presentan algunas recomendaciones atinentes al cuido del bienestar de los hogares.

Extending lagrange interpolation to develop a 4-node quadrilateral element

Today finite element method is a well established tool in engineering analysis and design. Though there axe many two and three dimensional finite elements available, it is rare that a single element performs satisfactorily in majority of practical problems. The present work deals with the development of 4-node quadrilateral element using extended Lagrange interpolation functions. The classical univariate Lagrange interpolation is well developed for 1-D and is used for obtaining shape functions. We propose a new approach to extend the Lagrange interpolation to several variables. When variables axe more than one the method also gives the set of feasible bubble functions. We use the two to generate shape function for the 4-node arbitrary quadrilateral. It will require the incorporation of the condition of rigid body motion, constant strain and Navier equation by imposing necessary constraints. The procedure obviates the need for isoparametric transformation since interpolation functions are generated for arbitrary quadrilateral shapes. While generating the element stiffness matrix, integration can be carried out to the accuracy desired by dividing the quadrilateral into triangles. To validate the performance of the element which we call EXLQUAD4, we conduct several pathological tests available in the literature. EXLQUAD4 predicts both stresses and displacements accurately at every point in the element in all the constant stress fields. In tests involving higher order stress fields the element is assured to converge in the limit of discretisation. A method thus becomes available to generate shape functions directly for arbitrary quadrilateral. The method is applicable also for hexahedra. The approach should find use for development of finite elements for use with other field equations also.